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The extragalactic background light (EBL) is all the electromagnetic energy released by resolved 
and unresolved extragalactic sources since the recombination era. Its intensity and spectral shape 
provide information about the evolution of galaxies throughout cosmic history. Since direct ob¬ 
servations of the EBL are very difficult to perform, the study of the interaction between the low 
energy EBL photons and high energy photons from distant extragalactic sources becomes rele¬ 
vant to constrain the EBL intensity. The main goal of this study is to investigate the opacity of 
the EBL to gamma rays by observing a sample of active galaxies nuclei (AGN) with the High Al¬ 
titude Water Cherenkov (HAWC) Gamma-Ray Observatory. Current gamma-ray observations up 
to 20 TeV performed by Imaging Atmospheric Cherenkov Telescopes (IACTs) have constrained 
the EBL intensity in the 0.1-50 pm region. HAWC, which monitors the gamma-ray sky in the 
100 GeV to 100 TeV energy range, will be able to detect about 10 AGN with the first year of 
HAWC data based on the extrapolation of steady-state spectra from the GeV band to TeV, and 
thus constrain the EBL in the poorly-measured 1-100 pm range. 
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1. Introduction 

The extragalactic background light (EBL) contains all the electromagnetic radiation released 
from resolved and unresolved galaxies since the recombination era. The EBL is the second most 
dominant radiation field in the universe, surpassed by the cosmic microwave background (CMB). 
Direct measurements of the EBL are difficult at most wavelengths, specially at ~ 15 jum where the 
contribution of the zodiacal light is some orders of magnitude larger than the EBL. Its spectrum 
is bimodal with one component peaking at ~ 1 /./m. which corresponds to stellar radiation, and a 
second component peaking at ~ 100 /./ m. created by the absorption of UV/optical light that is re¬ 
radiated by dust at infrared wavelengths. The EBL intensity and spectral shape, therefore, contains 
important information about the cosmic evolution of these astrophysical objects, the obscuring dust, 
and the relative contributions of starburst galaxies and AGN to the energy released over cosmic time 
[ 1 ]. 

The High Altitude Water Cherenkov (HAWC) Observatory is a ground-based TeV gamma-ray 
observatory in the state of Puebla, Mexico at an altitude of 4100 m a.s.l. Completed in March 
2015, the detector continuously measures the arrival time and direction of gamma-ray primaries 
within its 2 sr field of view. It is most sensitive to gamma-ray energies ranging from 100 GeV - 
100 TeV. Basic cuts on the data can be used to differentiate gamma-ray air showers from the large 
cosmic-ray background and employs cuts to differentiate them from cosmic-rays [2, 3, 4]. 

The main goal of this work is to select a sample of AGN to study the opacity of the EBL to gamma 
rays by observing in the 100 GeV to 100 TeV energy range. The spectra observed with HAWC will 
be corrected using an EBL model in order to estimate the the de-absorbed spectrum. By comparing 
the de-absorbed fluxes to the intrinsic spectra predicted by AGN models, we can constrain the EBL 
in the poorly-measured 1-100 jum region. 

2. Opacity to gamma rays 

On route to Earth, very high-energy gamma rays from cosmological sources have to pass 
through the radiation field of the EBL, resulting in their attenuation by pair producing interactions 
with photons in the near- to mid-IR regime of the EBL [5]. For a photon emitted by a source located 
at redshift z s and observed at z = 0 with observed energy Ey, the optical depth between z s and zo 

(0 < zo < Zs) is: 

tyy{E Y , zo,z s ) = J_ dz^ j j d/J. ~y~ d ' h ( f £ ’^ g yy( £ y( 1 + z) , e, 6 ), (2.1) 

where /./ = cos A, is the EBL proper photon density, ‘E = c\jz\, with l the proper distance, 6 
is the interaction angle between the two photons, Oyy{Ey{ 1 + z),£, 6) is the cross section of this 
interaction and £ f /, is the pair production threshold energy: 

2 mlc 4 

£th = ~r • (2-2) 

Ey{\ +Z)(1 - ft) 

The optical depth is a function of the observed energy and the source redshift, therefore the de¬ 
tection of extragalactic sources with HAWC is limited by redshift. If one expects to measure the 
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Figure 1: Relation between the observed energy and redshift for T = 1 (upper panel) and T = 3 (lower panel) 
predicted for 3 different EBL models ([7],[6],[13]) 


energy spectrum of a source at ~ 50 TeV then its redshift must be below 0.1 (Figure 1). The 
observed spectrum will be ~ 37% of intrinsic spectrum for T = 1 (equation 2.3). 


dN 

~dE 
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(2-3) 


The predicted opacity values given by [ 6 ] were interpolated in redshift and energy to find the 
redshift limit for HAWC’s observations and it was found that sources with z < 0.5 are mostly likely 
to be detected at very high energies with HAWC. 

For gamma rays of energy E-, the cross section for the 77 interaction peaks for low-energy photons 
with energy: 


2„4 
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where m e is the electron mass and c is the speed of light. In terms of wavelength [7]: 


(2.4) 


A [Atm] = 1.24£ 7 [TeV]. 

Since HAWC is sensitive between 100 GeV and 100 TeV, observations of AGN will allow us to 
study the EBL between 0.1 and 100 Atm. This is the range that corresponds mainly to the second 


3 





























AGN sample to study the EBL with HAWC 


Sara Coutino De Leon 1 


peak in the EBL and is dominated by IR emission from starlight absorbed and re-radiated by dust. 


3. AGN Sample 

To estimate the sensitivity of HAWC to AGN in the Northern Hemisphere, we selected AGN 
from the 1FHL catalog, the First Fermi-LAT Catalog of Sources above 10 GeV [8]. The 1FHL 
catalog contains 261 AGN in the field of view of HAWC, which decreases to 171 objects when we 
apply a redshift cut of z < 0.5. 

The Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope detects photons 
between 20 MeV and 300 GeV. Since the attenuation due to EBL is almost negligible in this energy 
range, we assume the energy spectra observed with the LAT are close to the intrinsic spectra of the 
sources [9]. We extrapolate the fluxes observed at GeV into the TeV band, and then attenuate the 
spectra using the EBL model of Franceschini et al. [6], Using the extrapolated and attenuated 
energy spectra, we can estimate the expected sensitivity of HAWC to the 171 1FHL objects at 
redshifts below 0.5. 

Based on the expected sensitivity of the full HAWC detector [2], we can estimate the time needed 
to detect the extrapolated AGN fluxes at the 5(7 level. Figure 2 shows the predicted TeV spectra 
of four AGN when compared to the 1-year sensitivity curve of the HAWC detector. Using the 
estimated one-year sensitivity of HAWC and simulated gamma rays from these sources, we found 
that 12 out of 171 AGN have a high probability of being observed with one year of data. These 
12 sources are listed in Table 1. From this list, six sources have variability periods registered with 
other telescopes. Measuring a spectra during one of the high flux states will provide more statistics 
in an EBL future work. 


3.1 Markarian 421 

Markarian 421 (Mrk 421) is one of the most studied gamma-ray AGN as well as the one 
that shows high degree of variability. It has been observed by HAWC, which found a possible 
correlation between the TeV and X-ray activity [10]. The fact that its distance from us is relatively 
short provides the opportunity to study the EBL using very high energy photons coming from it. 
As described in [3], HAWC data is divided in 10 bins, in each of them the flux normalization was 
calculated individually using the Likelihood Fitting Framework (LiFF), described in detail in [11]. 
It was found that the normalization flux depends very weakly on the assumed spectral index. Taking 
the median energy in each bin, preliminary energy flux measurements for a cut-off power law were 
calculated with 180 days of data during the construction phase of HAWC (2013/08 - 2014/03). 
Figure 3 shows that this preliminary flux measurement is in agreement with MAGIC observations 
[12]. Since there are very large bin-to-bin correlations, it is not a true differential energy spectrum, 
however we expect to obtain spectra with more accurate energy proxies to remove the correlations 
between bins. 
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Figure 2: Four of 171 1FHL extrapolated spectra to 100 TeV (red dots), the predicted attenuated spectra 
using and EBL model from [6] (blue dots), and the sensitivity curve for 1 year of observations with HAWC 
(black line) labeled in the upper left corner. The extrapolation is based on the steady-state spectra from the 
GeV band to TeV, and we have ignored the possibility of strong flares. 


4. Conclusion 


The selected extragalactic sources shown in Table 1 are the most likely to be detected within 
the first year of observations with HAWC, all of them have a redshift below 0.3, which according 
to the Figure 1, will be observed at E > 10 TeV. And over 180 days of data taken between August 
2013 and March 2014 with an early configuration of the HAWC Observatory have been analysed 
with a likelihood technique to produce flux measurements of Mrk 421. More days of data are 
needed to start an EBL formal study but with data taken with only a third of the detector array 
we have demonstrated the capability to detect Mrk 421 and the opportunity to study the EBL in a 
wavelength region that has not been constrained yet. 
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Figure 3: Preliminary energy flux measurements of Markarian 421 (black dots) for a cut-off power law in 
comparison with MAGIC observations [12] (red stars). The energy bins are the median energy from the 
energy distribution in each data bin where the normalization flux was fitted individually in each bin using 
the likelihood technique. Since there are very large bin-to-bin correlations, it is not a true differential energy 
spectrum, however we expect to remove the correlations between bins with a more accurate energy proxies. 
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1FHL 

Alternative name 

Type 

z 

r 

o/y/yr 

J0035.9+5950 

1ES 0033+595 

bzb 

0.086 

1.744=0.18 

6.02 

J0152.6+0148 

PMN JO 15 2+0146 

bzb 

0.080 

1.77 ±0.34 

4.85 

J0316.6+4119 

IC 310 

rdg 

0.019 

1.31 ±0.45 

13.16 

J0521.7+2113 

VER J0521+211 

bzb 

0.118 

1.97 ±0.14 

3.02 

J0650.8+2504 

1ES 0647+250 

bzb 

0.203 

1.56 ± 0.18 

10.25 

J0816.3-1310 

PMN J0816-1311 

bzb 

0.046 

2.06 ±0.27 

3.19 

J1104.4+3812 

Mrk 421 

bzb 

0.031 

1.91 ±0.06 

6.23 

J1230.8+1224 

M 87 

rdg 

0.004 

1.25 ±0.50 

20 

J1653.9+3945 

Mrk501 

bzb 

0.034 

1.86 ±0.10 

5.30 

J1728.3+5014 

IZw 187 

bzb 

0.055 

1.67 ±0.34 

3.85 

J2322.5+3436 

TXS 2320+343 

bzb 

0.098 

1.51 ±0.32 

9.68 

J2347.0+5142 

1ES 2344+514 

bzb 

0.044 

1.48 ±0.18 

5.14 


Table 1: Selected sources to study the EBL. The first column is the name in the 1FHL catalog, the second 
one is the alternative name given in the same catalog; the third column is the type of galaxy (bzb: blazar, 
rdg: radio galaxy), the fourth column is the redshift of the source, the fifth column is the spectral index in 
the 10-500 GeV regime given in the catalog, and the last column is the expected significance level of the 
detection during the first year of observations. 
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